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Abstract 

We develop a multifractal random tilling that fills the square. The multifractal is formed by 
an arrangement of rectangular blocks of different sizes, areas and number of neighbors. The 
overall feature of the tilling is an heterogeneous and anisotropic random self-affine object. The 
multifractal is constructed by an algorithm that makes successive sections of the square. At 
each n-step there is a random choice of a parameter pi related to the section ratio. For the case 
of random choice between pi and p2 we find analytically the full spectrum of fractal dimensions. 
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I - INTRODUCTION 

In the reference a deterministic multifractal tilling is developed to study percola- 
tion in multifractal systems. This object is a natural generalization of the square lattice 
in which it is possible to estimate analytically the full spectrum of fractal dimensions. 
Besides it exhibits a rich distribution of area among the blocks and a non-trivial topology, 
the number of neighbors varies along the object. The multifractal, however, is determin- 
istic in its construction. In spite of its properties it lacks the random ingredients found in 
natural systems. In order to make the model more realistic we develop a random version 
of this model. 

In this work we develop a multifractal that has the following characteristics. It has 
an infinite number of subsets each one with a distinct fractal dimension. For a specific 
case we determine analytically the dimension of each subset. The sum of all the subsets 
fills completely the square. The algorithm of construction of the multifractal has two or 
more free parameters, p^, that are related to the stretching of its subsets. During the 
construction of the multifractal there is random choice of pi which makes the object non 
deterministic. For the special case of just one parameter p = 1 the object degenerates 
into the square lattice. The extension of the model for an arbitrary dimension is trivial. 
Finally the algorithm for construction of the multifractal is simple and easily implemented 
in the computer. 

Geologic structures can show multifractal characteristics 

mm 

in physical 

properties like porosity, permeability, or sound velocity. In addition the rocks and strata 
are not deterministic objects, but present random features. These geologic systems that 
combine multifractality and randomness are natural candidates to be modeled by the 
proposed model. 

The deterministic multifractal is an intuitive generalization of the square lattice. Sup- 
pose that in the construction of the square lattice we use the following algorithm: take a 
square of size 1 and cut it symmetrically with vertical and horizontal lines. Repeat this 
process /-times; at the l^^ step we have a regular lattice with 2^ x 2^ cells. The setup algo- 
rithm of the deterministic multifractal is quite similar, the main diflFerence is that we do 
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not cut the square in a symmetric way, but according to a given constant ratio. The ran- 
dom multifractal is an extension of the deterministic model. In the random multifractal 
we use a non constant ratio in the algorithm setup. 

The main objective of this work is to present a random multifractal tilling and de- 
termine its set of fractal dimensions for a specific case. The manuscript is organized as 
follows: in section 2 we present the deterministic multifractal and review the derivation 
of its set of fractal dimensions. In section 3 we develop the random multifractal and in- 
dicate a formula that estimate its set of fractal dimensions for the case of random choice 
between two given distinct pi. Finally in section 4 we summarize the main diflFerences 
between deterministic and random multifractal tillings and give our final remarks. 



II - THE DETERMINISTIC MULTIFRACTAL 

We start with a square of linear size 1 and a parameter < p = ^ < 1, for integers 
r and s. The first step, n = 1, consists of making one vertical section of the square. 
The square is cut in two pieces of areas ^ = and ^ = by a vertical segment. 
This process is shown in figure (a), where we use p = |. The step n = 2 consists of a 
horizontal cut with the same section p, this step is shown in figure (b). At this point 
there are four rectangular blocks: the largest one of area (j^)^? two of area (j^p- and 
the smallest one of area (j^)^- We call a set of all elements with the same area as a 
/c-set. At the step n = 2 there are three fc-sets. 

In the third step, n = 3, we repeat the same process of vertical section as in step n = 1. 
At n = 4 we make new horizontal sections. All the sections have the same ratio p. At 
the n^^-step there are 2^ blocks. The partition process produces a set of blocks with a 
variety of areas. At the step n we have n + 1 fc-sets. The construction process of the 
deterministic multifractal is done by vertical and horizontal sections of the square. As the 
section ratio is always the same and the segment section is displaced in space the object 
shows self- affinity. 

At the n^^-step of the algorithm, the partition of the area of the square follows the 



3 



binomial rule: 



E c*," ( (^^j - (^^j , (1) 



or: 

The number of elements of each /c-set in this case is C^. 

Each fc-set determines a specific dimension. For the estimation of the spectrum of 
fractal dimension Dj. of an object X we use the box counting method [8]. The object X is 
immersed in the plane of real numbers 3?^, we use the trivial metric. We cover the object 
by just-touching square boxes of side length e. Let N{X) denote the number of square 
cells of side length e which intersect X. If 

log N{X) 

Dx = hm.^Q — — , (3) 

is finite, then Dx is the dimension of the X. 

In our case the object X is a fc-set. Remember that the fc-set corresponds to a set of 
rectangles of the same area. If the square is partitioned n times we have e = l/{s + r)^. 
For each k-set the area of blocks (using e area units) is done by: 

Nk = C^ 5^r("-^\ (4) 

where is the binomial coefficient that express the number of elements fc-type, and 
rj^{n-k) -g ^Yie area of each element of this set. We put together all these elements to 
have the fractal dimension of each fc-set: 

D,-Um J'^^'-Hm "^""'^ (5) 

log^ log [s + r)^ 

The case r = 5 = 1 is degenerated. In this situation the subsets of the lattice are composed 
uniquely by square cells of the same area. Therefore the tilling is formed by a single subset 
of dimension 2. 

Figure [21 shows the spectrum of for several values of p. We plot in the x-axis the 
normalized fc-set, fc/n, and in the y-axis the respective D^. The dashed line corresponds 
to a p that is very closed to 1 in this situation the blocks have almost the same area and 
Dk is a symmetric curve. In the opposite situation p ^ we have curves that correspond 



to a tilling in which the blocks are very stretched in one direction. In the figure it is also 
depicted the curves corresponding to p = 0.4, 0.6 and 0.8. 

Ill - THE RANDOM MULTIFRACTAL 

In the last section we use a constant p through all the n-steps of the construction of 
the multifractal. A more realistic generalization of the above process is to use diflFerent 
values of the ratio p along the growing process of the multifractal. In the general case a 
random pi is used at each step of the algorithm. In this work we explore in more detail 
the case where at each n-step of the construction of the multifractal we randomly choose 
between two rational values pi = Si/ri and p2 = S2/r2 for ri + Si = r2 + S2. The picture 
of this process is similar to figured but this time the section of the square does not follow 
a constant ratio making the partition of the square asymmetric. 

We have in the algorithm of the random multifractal the following choices. At n = 1 
we divide the square in two blocks. There are two partition possibilities for the square: 
pi/(l+pi) and (l-pi)/(l+pi) or ^2/(1+^2) e (l-p2)/(l+P2)- At the second step, n = 2, 
there are two possibilities for each former configuration. These are the three partition 
possibilities at step n = 2: 

+ Pl)^ - + (1 - Pi)V(i + Pi)', Pi(i - Pi)/(i + Pif: 

or piP2/(l + Pi)(l + P2), (1 - P1)P2/(1 + Pi)(l + P2), Pi(l - P2)/(l + Pi)(l + P2), 
(1-Pl)(l-P2)/(1+Pl)(l + P2), 

or pl/il + p,)\ p,{l - P2)/(1 + P2)^ (1 - P2)V(1 + P2)^ P2(l - P2)/(l + 92?. 

Note that there are not four configurations but only three. In fact, it is the same to 
make a section pi after p2, or P2 after pi. 

At the step n, pi is chosen / times and p2 is chosen m times for / + m = n. At each 
step the areas of the blocks are obtained by the multiplication of the area of the block of 
the former step. Therefore the partition of the area of the square is done by the following 
equation: 

■=55<rf.)'(IS)'>(rfty(feSy" 
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This formula is analogous to equation ^ of the deterministic case. The occurrence of 
blocks of the same fc-sets (objects of same area) is smaller for this case compared with the 
deterministic model. Otherwise, the quantity of fc-sets is larger, it is done by (/ + l)(m+l). 

The determination of the spectrum of fractal dimensions for the case of the multifractal 
is similar to the deterministic case, equation The equation for in this case is the 
following: 

log Cls\r[-' C';'sir^-' 

Dk = hrun^oo — 7 , (7) 

log {si + riY 

where the set of k corresponds to any combination of / and m such that n = I + m. The 
generalization of this formula for several pi (for the case where each pi obeys ri + Si = 
Ti + Si) is performed by the introduction of new terms CJ^' S2r^'~^ in the numerator. 

Figure El shows the curve in crescent order of fractal dimensions D^- In the ^-axis is 
plotted Dk and in the x-axis, {k/nY ^ the respective fc-set. We notice that the /c-sets in 
the X-axis are shown in order of crescent D^. The area of the fc-sets is normalized between 
zero and one. In the figure there are three curves pi = 1/4 (point line), p2 = 2/3 (dashed 
line), and the random curve (thick line). In this simulation we use n = 100. The random 
curve is defined as the spectrum that results from the random choice between pi and p2- 
The estimation of the random curve is performed assuming that for n ^ oc we have 
m = j = n/2. 

A typical view of the the deterministic and the random multifractals is shown in figure 
H These two pictures show the deterministic multifractal (a) and one realization of the 
random multifractal (b). In figure 0] (a) we use p = 1/5 and in figure 0] (b) pi = 1/4 and 
P2 = 2/3. In both situations we evolved the algorithm until n = 10. 



IV - FINAL REMARKS 

In this work we develop a random self-affine multifractal object that forms a tilling of 
the square. The multifractal is constructed by an algorithm that make successive sections 
of the square. In the special case where the section ratio is 1/2 the algorithm generates the 
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well known square lattice. At each n-step there is a random choice of partition parameter 
Pi. We explore in detail the case in which the ratio choice done is between two specific 
values pi and p2- For this situation at the limit of infinite partition of the square, n ^ oc, 
we find analytically the full spectrum of fractal dimensions. 

In the references the deterministic multifractal object was used to study 

percolation properties and the breaking of universality class in two dimensions. In fact, 
the critical exponent related to the correlation length for the deterministic multifractal is 
not universal, it means, it is not the same as in standard percolation in two dimensions. 
This fact is related with the breaking of local isotropy. There are some points to consider 
in the symmetry properties of this multifractal object and its isotropy. The first is related 
with the stretching of the blocks of the multifractal which increases as p ^ 0. Moreover 
the topology of multifractal, specially its coordination number, changes along the object. 
The object we develop in this work is a natural candidate to extend the cited papers to 
a non deterministic analyses of the breaking of universality class. In special we intend to 
analyze the reason for the breaking of universality class in the random multifractal and 
how the local isotropy affects this process. 



The authors gratefully acknowledge the financial support of Conselho Nacional de 
Desenvolvimento Cientifico e Tecnologico (CNPq)-Brazil, FINEP and CTPETRO. 



[1] G. Corso, J. E. Freitas, L. S. Lucena, R. F. Scares, accepted in Phy. Rev. E (2003), 
|cond-mat/0212530l 

[2] L. S. Lucena, J. E. Freitas, G. Corso, and R. F. Scares, Brazilian Journal of Physics 33 637 
(2003). 

[3] J. Muller, Annales Geophysicae - Atmospheres Hydrospheres and Space Sciences 11 (6): 525 
(1993). 

7 



[4] P. N. Khue, O. Fluseby, A. Saucier, J. Muller, Journal of Physics - Condensed Matter 14 
(9): 2347 (2002). 

[5] R. H. Riedi, Multifractals and wavelets: a potential tool in Geophysics, Proceedings of the 

68th SEG Meeting, New Orleans, Louisiana, USA, (1998). 
[6] P. Hubert, Hydrological Sciences Journal 46 897 (2001). 
[7] S. Sachs, S. Lovejoy, D. Schertzer, Fractals 10 253 (2002). 
[8] M. Barnsley, Fractals Everywhere, (Academic Press, Boston, 1988). 
[9] J. E. Freitas, G. Corso, and L. S. Lucena, |coird-mat/03107"79l (2003). 

FIGURE LEGENDS 

FIG. 1: The initial two steps in the formation of the deterministic multifractal. In (a) a vertical 
line cut the square in two pieces according to p. Two horizontal lines sectioning the rectangles 
by the same ratio are depicted in (b). The same process is repeated n-times inside each block. 
At the limit n ^ oc we have a multifractal. 



FIG. 2: The spectrum of fractal dimensions, Dj^, for the deterministic multifractal. In the y-axis 
is Dj^ and in the x-axis the respective normalized fc-set. The curves for several p are indicated 
in the figure. 



FIG. 3: The spectrum Dj^ for pi = 1/4 (point line) and p2 — 2/3 (dashed line). In thick line is 
shown the random curve. In the the x-axis we show {k/nY the fc-sets in order of crescent D]^. 
The area of the fc-sets are normalized to one. 
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FIG. 4: The random multifractal tilling evolved until n = 10 for (a) the deterministic multifrac- 
tal (p = 1/3) and (b) one realization of the random multifractal (for pi = 1/4 and p2 = 2/3). 
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